﻿using System;
using System.Text;
using System.Drawing;
using System.Buffers;
using System.Collections;
using System.Collections.Generic;
using System.Runtime.InteropServices;

public static partial class glDRIVE
{
    /*
    函数 gl.simp
    Simpson求积法
    f计算被积函数值f(x)的函数名。
    参数 a: 积分下限。
    参数 b: 积分上限。要求b>a。
    参数 eps: 积分精度要求
    返回值 函数返回积分值。
    */

    public static string drive_simp()
    {
        double a, b, eps, t;
        a = 0.0;
        b = 1.0;
        eps = 0.000001;
        t = gl.simp(a, b, eps);

        return gl.html_table("Simpson积分[y=Math.Log(1.0 + x) / (1.0 + x * x)]：", new double[3] { a, b, t });
    }

    // 计算被积函数值
    private static double simpf(double x)
    {
        double y;
        y = Math.Log(1.0 + x) / (1.0 + x * x);
        return (y);
    }
}